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Abstract 

Neutrinoless double-beta decay nuclear transition matrix elements are generated by an effective two-body 
transition operator and it consists of Gamow-Teller like and Fermi like (also tensor) operators. Spectral 
distribution method for the corresponding transition strengths (squares of the transition matrix elements) 
involves convolution of the transition strength density generated by the non-interacting particle part of the 
Hamiltonian with a spreading function generated by the two-body part of the Hamiltonian. Extending 
the binary correlation theory for spinless embedded /c-body ensembles to ensembles with proton-neutron 
degrees of freedom, we establish that the spreading function is a bivariate Gaussian for transition operators 
0{ko) that change ko number of neutrons to ko number of protons. Towards this end, we have derived 
the formulas for the fourth-order cumulants of the spreading function and calculated their values for some 
heavy nuclei; they are found to vary from ~ —0.4 to —0.1. Also for nuclei from 76 Ge to 238 U, the bivariate 
correlation coefficient is found to vary from ~ 0.6 — 0.8 and these values can be used as a starting point for 
calculating nuclear transition matrix elements using the spectral distribution method. 
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I. INTRODUCTION 



Double-/? decay (DBD) is an extremely rare weak-interaction process in which two identical 
nucleons inside the nucleus undergo decay with or without emission of neutrinos. Theoretically, 
the two neutrino double beta decay (2vf3~f3~) process was first predicted by Mayer [1] following 
the suggestion of Wigner. This process has been observed in more than 10 nuclei and best adopted 
values for the half-lives were tabulated recently by Barabash [2] . In 1937, following the suggestion 
of E. Majorana [3], Racah [4] pointed out the possibility of neutrinoless double-/? decay (NDBD 
or 0v/3~/3~). Furry [5] in 1939, for the first time calculated NDBD half-lives. Fundamental sig- 
nificance of NDBD is that its experimental confirmation will tell us about lepton number violation 
in nature and that neutrino is a Majorana particle. More importantly, NDBD gives a value or a 
bound on neutrino mass [6] provided the half-lives are known experimentally and the correspond- 
ing nuclear transition matrix elements (NTME) are obtained using a reliable nuclear model. So 
far only Klapdor et al [7] claimed to have evidence (at a confidence level of 4.2<r) for 0is/3~/3~ in 
76 Ge. At present large number of NDBD experiments are being carried out and many others are in 
development stage in various laboratories around the world. The nuclei being considered are 48 Ca, 
76 Ge, 82 Se, 100 Mo, 116 Cd, 130 Te, 136 Xe, 150 Nd and so on [6]. Following this, several nuclear models 
are employed for calculating the NTME for various candidate nuclei mentioned above. Some of 
the models used for NDBD studies are shell model using recent state-of-the-art large scale calcu- 
lations [8, 9], quasi-particle random-phase approximation with various extensions [10, 11], inter- 
acting boson model [12], pseudo-SU(3) model [13], projected Hartree-Fock-Bogoliubov method 
including pairing plus quadrupole-quadrupole interaction [14, 15], generating coordinate method 
with particle number and angular momentum projection [16]. 

Statistical spectral distribution theory [17] gives a method for calculating transition strengths 
(squares of transition matrix elements) generated by a transition operator. This theory starts with 
shell model spectroscopic spaces and the same shell model inputs (single particle energies and 
effective two-body interactions). Here one constructs smoothed forms (spectral distributions) for 
various observables ignoring the fluctuations and this is based on random matrix representation of 
Hamiltonians (also other operators), unitary decompositions of the operators and quantum chaos. 
Spectral distribution theory has been applied in the past, with various approximations, to a variety 
of problems in nuclear structure and they include (i) bound on time-reversal non-invariant part 
of the nucleon-nucleon interaction [18], (ii) single particle transfer [19], (iii) /3-decay rates for 
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pre-supernovae evolution [20, 21], (iv) giant dipole widths [22], and (v) parity breaking matrix 
elements in compound nuclei [23]. Our purpose in this paper is to describe spectral distribution 
method for the NTME for NDBD and establish that the spreading function that enters in the theory 
is close to a bivariate Gaussian. With this result, it is possible to apply in future spectral distribution 
method to NDBD. Now we will give a preview. 

In Sec. H, we describe briefly the structure of the Oz//3~/3~ transition operator which is a two- 
body operator. Section III describes the spectral distribution method for NDBD. Section IV gives 
the binary correlation theory for traces of operators over two-orbit (proton and neutron) configu- 
rations. Using these, in Sec. V, derived are the results for the fourth-order bivariate cumulants and 
the bivariate correlation coefficient for the spreading transition strength density function appro- 
priate for NDBD. Numerical results are presented for the fourth-order cumulants to demonstrate 
that the transition strength density is close to a bivariate Gaussian. In addition, for several heavy 
nuclei, the values for the bivariate correlation coefficient are presented. Finally, Sec. VI gives 
conclusions. 



II. 0u DOUBLE-BETA DECAY TRANSITION OPERATOR 

Half-life for 0u double-beta decay, for the 0+ ground state (gs) of a initial even-even nucleus 
decay to the 0^ gs of the final even-even nucleus is given by [24] 

[nm opr 1 = <?* im^i 2 f^V , (d 

V m e J 

where (ra„) is the effective neutrino mass (a combination of neutrino mass eigenvalues and it also 
involves neutrino mixing matrix). The G 0u is phase space integral (kinematical factor) dependent 
on charge, mass and available energy for NDBD process; tabulations for G 0u are available [25, 26]. 
The M 0u is the NTME generated by the NDBD transition operator and it is a sum of a Gamow- 
Teller like (M GT ), Fermi like (M F ) and tensor (M T ) two-body operators. As it is well known that 
the tensor part contributes only up to 10% of the matrix elements [27, 28], we will neglect the 
tensor part. Then we have, from the closure approximation which is well justified for NDBD [24], 



M 0u (0 + ) = M^(0+) - 4 m f (0 + ) = (0/ II 0{2 : Ou) || 0+) , 

9a 

O(2:0i/) = Y. n ^) T a<{ a ^^- g T) ■ 

a,b V 9aJ 
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As seen from Eq. (2), NDBD half-lives are generated by the two-body transition operator 0(2 : 
Ov); note that a, b label nucleons. The qa and g v are the weak axial-vector and vector coupling 
constants. The "H(r aft , E) in Eq. (2) is called 'neutrino potential'. Here E is the average energy of 
the virtual intermediate states used in the closure approximation. The form given by Eq. (2) is jus- 
tified only if the exchange of the light Majorana neutrino is indeed the mechanism responsible for 
the NDBD. The neutrino potential, defining completely the two-body NDBD transition operator 
0(2 : Ou) is, to a good approximation, given by [27-31], 

U(r ab ,E) = — ®(r ab ,E). (3) 

Tab 

Here, R in fm units is the nuclear radius and similarly r ab is in fm units. The function $ is given 
by [27, 28,31], 

$(r a6 ,£) = - 

7T 

In Eq. (4), fi(x) = — f^t -1 cost dt = Ci(x) = 7 + In x + J* t" 1 (cost — 1) dt and fi{x) = 
— j™t~ l s\ntdt = Si(x) — |; Si(x) and Ci(x) are the sine and cosine integrals. It is useful 
to mention that $(r ab ,E) ~ exp(-|£r afe ). Note that Tic = 197.327 MeV fm. The effects 
of short-range correlations in the wavefunctions are usually taken into account by multiplying 
the wavefunction by the Jastrow function [1 — 7 3 e~ 7ir «i>(l — 72^)] [28, 29, 31]. Now keeping 
the wavefunctions unaltered, the Jastrow function can be incorporated into TL(r ab: E) giving an 
effective n eff (r ab ,E), 

U(r ab ,E) -+ n eff (r ab ,E) = H(r ab ,E)[l - 73 e~ 71 r -»(l - 72 r 2 ab )} 2 . (5) 

The standard set of values for the parameters 7 l5 72 and 73 are given ahead. The most important 
point about Eq. (1), as mentioned earlier, is that an experimental value of (bound on) T^ 2 will 
give a value for (bound on) neutrino mass via Eq. (1) provided we know the value of the square 
of NTME M 0u {f) + ) of the NDBD two-body transition operator 0(2 : Ou), connecting the ground 
state of the initial and final even-even nuclei involved. 

Let us say that for the nuclei under consideration, protons are in the single particle (sp) or- 
bits j p and similarly neutrons in j n . Using the usual assumption that the radial part of the 
sp states are those of the harmonic oscillator, the proton sp states are completely specified by 
(n p ,l p , j p ) with ii p denoting oscillator radial quantum number so that for a oscillator shell J\f p , 
2n p +£ p = Af p . Similarly the neutron sp states are (n n , £ n , j n ). In terms of the creation (a 1 ") and an- 
nihilation (a) operators, normalized two-particle (antisymmetrized) creation operator Aj l (jij 2 ) = 
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)J and then |0) = \(jij 2 )J^) is a normalized two-particle state. With 
these the NDBD transition operator can be written as, 



0{2 :0u)= °b*J5,?M 5>Z(J?#) K(^?)} f 

fi>3%;jS>32;J 1 1 



(6) 



Note that £>i . P . .„ .„(0i/) = {{j\il)JM | 0(2 : Ou) \ {j%jl)JM) a are two-body matrix elements 
(TBME) and 'a'denotes that we are considering antisymmetrized two-particle wavefunctions; J 
is even for j 1 = j 2 or j 3 = j 4 . Numerical values for the TBME 's in Eq. (6) follow from the 
definition of 0(2 : Ou) in Eq. (2) which can be expressed as 



0{2 : Ou) = Ea, b r+r+(0 GT - {g 2 v /g A )0 F ) ; 
O gt = a a ■ a h H eff (r ab , E) , O f = H eff (r ab , E) . 



(7) 



Note that t+t 6 + simply changes two neutrons to two protons and for the remaining parts, to obtain 
the TBME, we use Brody-Moshinsky brackets [32-34]; see [27] for an alternative approach. Then, 
the TBME are given by, 
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Here, } are the 9j-coefficients, (• • • | • • •) are Brody-Moshinsky brackets [32-34], B(- • • ) 

are Brody, Jacob and Moshinsky coefficients [33] and I p are Talmi integrals [34]. It is important to 
mention that antisymmetrized matrix elements for 0(2 : Ou) are used in the shell model calcula- 
tions of NTME while in QRPA related studies non-antisymmetrized matrix elements are employed 
[28, 35]. 
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There are a number of parameters in the NDBD transition operator and some of them are usu- 
ally varied in the calculations. The parameters are: (i) R = 1.1A 1 / 3 — 1.2A 1 / 3 fm [31]; (ii) 
b = 1.003A 1 / 6 fm [27]; (iii) E = 1.12A 1 / 2 MeV [30]; (iv) g A /g v = 1 (quenched) or 1.25 
(unquenched) [31]; (v) three different choices for the parameters (71, 72, 73) in Eq. (5) are 
(1.1,0.68, 1) [Miller-Spencer], (1.52, 1.88,0.46) [CD-Bonn] and (1.59, 1.45,0.92) [AV18]; these 
values are taken from [28, 31]. It is useful to mention that the kinematical factor G 0u depends on 
the coupling constant g& (standard value is 1.25) and also some calculations use different values 
for r in R = r^A 1 ^ fm; the standard value is r = 1.2. Then a scaling for M 0u is [31, 36] 



Now we will give the spectral distribution formulation for calculating NTME for NDBD. 

III. SPECTRAL DISTRIBUTION METHOD FOR NTME 

Following Eqs. (1) and (2) for NDBD half-lives, the corresponding NTME |M 0jy | 2 can be 
viewed as a transition strength (matrix element connecting a given initial state to a final state by a 
transition operator) generated by the two-body transition operator 0(2 : 0u). Therefore, spectral 
distribution theory [37-41], based on random matrix theory, for transition strength densities, can 
be applied [42]. Transition strength density is defined as the transition strength multiplied by the 
state densities at the initial and final energies involved. 

Let us consider shell model spherical orbits with angular momenta jf, jf, . . . with m p 
protons distributed over them and similarly m n neutrons over j™, . . . , j™ orbits. Then the 
proton configurations m p = [m p: m p: . . . ,m r p ] where m p is number of protons in the orbit jf 
and £[ =1 m p = m ~p- Similarly the neutron configurations m n = [m^, m^, . . . , to*] where 
m l n is number of neutrons in the orbit jf and ^=1 m n = m n- With these, (toJ,,to^)'s de- 
note proton-neutron configurations. The nuclear effective Hamiltonian is one plus two-body, 
H = h + V and we assume that the one-body part h includes the mean-field producing part 
of the two-body part. Thus [43, 44], V is the irreducible two-body part of H. The state den- 
sity I H (E), with (( )) denoting trace, can be written as a sum of partial densities defined 



over(m p ,TO n ),i.e. = ((S(H - E))) {m ^ > = £(^,550 ((6(H - E)))^^ = 

£(m p ,J^) I^'^KE) = £ (rrpi ^) d(m p ,m n )p^^(E). Here d(m p ,m n ) is dimension and 



p( m f" m ™)(£') is the partial density normalized to unity. For nuclear Hamiltonians, it is well under- 




(9) 
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stood [17, 40, 45] that the smoothed form for the eigenvalue densities is generated by the action 
(locally) of EGOE(l+2) [embedded GOE of one plus two-body interactions] operating in the 
Gaussian domain. This gives, 

j(m p ,m n )^ = £ lfi^{E) . (10) 

(m p ,m^) 

In Eq. (10), Q denotes Gaussian. The Gaussian partial densities are defined by the cen- 
troids E c (m p ,m n ) = (H)^^ ~ (h) ( ™^ and variances a 2 (m p ,m n ) = (H 2 )^'^ - 
[E c (m p , m n )} 2 ~ ^v 2 )' mp,m "' and, as they are traces over (m p , m n ) spaces, they can be calculated 
without recourse to H matrix construction. Propagation equations for them, in terms of the sp en- 
ergies defining h and TBME defining V, are easy to write down. Unitary group tensorial structure 
of h and V gives further simplifications of Eq. (10); see [43, 44, 46] for details and applications 
to /p-shell and also for heavy (A ^ 150) nuclear data analysis. 

Random matrix theory, based on EGOE(l+2), for the (smoothed) transition strength densi- 
ties I (Ei,Ef) = I(Ef) \ (Ef | O | Ei) \ 2 I{Ei) allows us to write I a as a convolution of the 
corresponding density generated by the mean-field part h with a spreading bivariate Gaussian 
Puv-q-.o-n due to the interaction V. This gives [37-40], for the square of the O matrix elements, 
with |0+) = \Ei) and |0^ } = \Ej) where E's are energies, 

\(E I O I E-)\ 2 = V Q — 9 - s ' 

{m p ,m n )i,(m p ,m n ) f 

(ID 

x \ {(m p ,m n ) f | O | (m^nQi)] 2 blv ' G: ° :V fj. ' c ' c ' - / ' — . 

(9j ''(fij) Pg ' f (Ef ) 

In Eq. (11), \((irf p ,m n )f \ O \ (m p ,m n )i)\ 2 is the mean square matrix element of O connecting 
(rn p , m n )i and (m p , m n ) / configurations, 

\{(m p ,m n )f | O | (m p ,m n )i)\ 2 = {d[(m p , m n )i\ d[(m p , nQf]}" 1 

x ^2 \((m p ,m n ) f /3 \ O \ (m p ,m n )i a}\ 2 . 

(12) 

For later reference it is also useful define, 

([(9]t(9)^^ = d[(m p ,m n ) f ] £ \((m p ,m n ) f | O \ (m P ,m n ) t )\ 2 . (13) 



Following [37-39], the two centroids E l c and E{ and the two variances of and a 2 of the marginal 
densities of the bivariate Gaussian pbiv-g-.o-.v, with some approximations are given by, 

E' l c = E c ((m p ,m n )i) , E{ = E c ((m p ,m n ) f ) , 

(14) 

of = o 2 ((m p ,m n )i) , aj = a 2 ((m p , m n ) f ) . 

These are the proton-neutron configuration centroids and variances generated by H. Eq. (14) has 
its basis in Eq. (54) given ahead. Although the general structure and importance of the bivariate 
correlation coefficient ( in Eq. (11) is well understood [37-39], an expression for ( in terms of 
traces over (m p , fn n )i and (m p , m n ) f spaces is not yet available. Then the standard approximation, 
for completing the theory, is 

, *n [K,m ra ) t ] ([0]tvQ0yP)^ m 

C= i =\ X PQ [{m pi m n ) t }= ^ ( m ;, m „ )i • (15) 

yXao [(m p ,m n )i\ X 02 [(m p , m n )i] \M u l 

To proceed further, propagation formulas for the traces in Eqs. (12) and (15) are needed. For the 
trace in Eq. (12), using the results in [50], it is easy to write the propagation formula in terms of 
the TBME O^_(0u) defined by Eq. (6), 

\{(m p ,m n ) f | O | (m p ,m n ),j)\ 2 d[(m p ,m n ) f ] 

v m^(q)K(/3) - 6 af) ][NM - - rnffi) - M 

' a ^ tS N n (a)[NM - 5^(^(5) - 5 lS ] 

(16) 

xE [^ w (0^)] 2 (2Jo + l), 

Jo 

(m p ,m n ) f = {m^rrQi x (l+ 1+1^1^) . 

Note that in Eq. (16) the final configuration is defined by removing one neutron from orbit a 
and another from (3 and then adding one proton in orbit 7 and another in orbit 5. Also, N p (a) 
is the degeneracy of the proton orbit a and similarly N n (a) for the neutron orbit a. In the limit 
Pbiv/pp = 1 in Eq. (11), substituting the result of Eq. (16) in Eq. (11) gives the NTME for 
NDBD. However, this zero-th order approximation will not be good as in general it is expected 
that ( > 0.5; see Sec. V ahead. Before proceeding to implement the theory given above, it is 
essential to test the important approximation used in the theory, i.e. the bivariate Gaussian form for 

8 



the spreading function generated by V by calculating the fourth order bivariate cumulants (they 
will be zero for a bivariate Gaussian). In addition we also need an expression for the bivariate 
correlation coefficient (. To address these two problems and provide generic results, we consider 
the spreading function defined over proton-neutron spaces i.e., 



where H is a two-body Hamiltonian. We will consider Hamiltonians that preserve (m p , m n ) and 
then H = H pp + H nn + H pn . This is quite appropriate for heavy nuclei. Adopting the binary 
correlation theory, the bivariate moments of p^^^ t ' ijnv,mn ^ f are evaluated by considering ran- 
dom /o-body H operators. Similarly, we represent the transition operator O by random /c^-body 
operator that changes k number of neutrons to k G number of neutrons. This is equivalent to 
using EGOE representation for both H and O operators [18, 40]. Let us mention that, from now 
onwards, we consider only the two-orbit configurations (m p , m n ) [for generality, these are denoted 
as (mi,m 2 ) in the next two sections]. 

IV. BINARY CORRELATION RESULTS FOR RANDOM HAMILTONIANS 

Binary correlation theory for moments defined over a single unitary orbit is given by Mon and 
French [47, 48] and they correspond to the moments generated by spinless EGOE(A;) in the dilute 
limit (dilute limit is defined in Sec. IV). The theory is extended to certain two-orbit moments by 
Tomsovic [49] . For the two problems mentioned in Sec. Ill, we need traces defined over two-orbits 
(protons and neutrons) with the H preserving the two-orbit symmetry and the transition operator 
O acting on a two-orbit configuration generating a unique final two-orbit configuration. In the 
present section, we will give the basic binary correlation results adopted for this situation and in 
Sec. V, we will consider their applications. 

A. Results for single unitary orbit 

Let us begin with a /c^-body operator, 



a, 

Here, a\k H ) is the k H particle creation operator and (5(k H ) is the k H particle annihilation op- 
erator. Similarly, v^f are matrix elements of the operator H in k H particle space i.e., v^f = 

9 



(m v ,m n )i, {m p ,m„)f 
Pbiv.O-.H 



(E i: E f ; E l c ,E f c ,a h a f ,C) , 




(17) 



{kn/3 | H | knOi) (it should be noted that Mon and French [47, 48] used operators with daggers 
to denote annihilation operators and operators without daggers to denote creation operators). Fol- 
lowing basic traces will be used throughout, 

/ \ m 

X)a t (*)a(*)=Q) ^r£a\k)a(k)\ = (™) . (18) 
/ \ m ~ 

^at(fc)S(A/)a(fc)= (^W) 

(20) 

^«(A;)5(A;V(A:)=( iV "^ A;/ )5(A: / ) 

(21) 



Equation (18) follows from the fact that the average should be zero for m < k and one for m = k 
and similarly, Eq. (19) follows from the same argument except that the particles are replaced by 
holes. Equation (20) follows first by writing the AZ-body operator B(k') in operator form using 
Eq. (17), i.e., 

B(k') = J2< ^(k'h(k'), (22) 
and then applying the commutation relations for the fermion creation and annihilation operators. 



This gives Y1b- 1 v b' Y^ a a\k)a(k)^(k'). Now applying Eq. (18) to the sum involving 



a. 



gives Eq. (20). Eq. (21) follows from the same arguments except one has to assume that B(k') is 
fully irreducible v — k' operator and therefore, it has particle-hole symmetry. For a general B(k') 
operator, this is valid only in the N — > oo limit. Therefore, this equation has to be applied with 
caution. 
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Using the definition of the H operator in Eq. (17), we have 



(H(k H )H(k H )) m = J2{ v h} (u\k H )P(k H )p\k H )a(k H )) T 



= v 2 H fe^N few'N «N (23) 



= v 2 H T(m, N, k H ) . 

Here, H is taken as EGOE(A;#) with all the kn particle matrix elements being Gaussian variables 
with zero center and same variance for off-diagonal matrix elements (twice for the diagonal matrix 



elements). This gives (v^f) 2 = v 2 H to be independent of a, (3 labels. It is important to note that 
in the dilute limit, the diagonal terms [a = f3 in Eq. (23)] in the averages are neglected (as they 
are smaller by at least one power of l/N) and the individual H's are unitarily irreducible. These 
assumptions are no longer valid for finite- N systems and hence, evaluation of averages is more 
complicated. In the dilute limit, we have 

T(m,N,k H ) = (xy(*ff) Ixi^m^imI^m 



m + kn\ ( m 
k H ) \k H 

Note that, we have used Eq. (19) to evaluate the summation over f3 and Eq. (18) to evaluate 
summation over a in Eq. (24). In the 'strict' N — > oo limit, we have 

In order to incorporate the finite- N corrections, we have to consider the contribution of the diago- 
nal terms. Then, we have, 



T(m, N, k H ) - 
Now we will turn to the fourth order averages 



m 
k H 



m + ku 
k H 



(26) 
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For averages involving product of four operators of the form (H(k H )G(k G )H(kH)G(k G )) m , 
with operators H and G independent and of body ranks k H and k G respectively, there are two 
possible ways of evaluating this trace. Either (a) first contract the H operators across the G oper- 
ator using Eq. (21) and then contract the G operators using Eq. (20), or (b) first contract the G 
operators across the H operator using Eq. (21) and then contract the H operators using Eq. (20). 
However, (a) and (b) give the same result only in the 'strict' iV — > oo limit and also for the result 
incorporating finite N corrections as discussed below. In general, the final result can be expressed 
as, 



H(k H )G(k G )H(k H )G(k G )} m = vjj v G F(m,N,k H ,k G ) . (27) 



In the 'strict' dilute limit, F(m ) N, kn, k G ) is given by 



k G J \k H J \k H J \k G/ 
In order to obtain finite-iV corrections to F(- • • ), we have to contract over operators whose lower 

symmetry parts can not be ignored. The operator H(kn) decomposes into irreducible symmetry 

parts denoted by s — 0, 1 , 2, . . . , kn with respect to the unitary group SU (N) . For a -body 

number conserving operator [48, 50], we have 

Here, the T{s) are orthogonal with respect to m-particle averages, i.e., (j r (s)J r ^(s')) m = 5 SS /. 
Now, (H(k H )G(k G )H(k H )G(k G )} m will have four parts, 



(H(k H )G(k G )H(k H )G(k G )) r 



a,/3,-y,6 

+v 2 H v G <« t (^)a(^)7 t (A; G )5(A; G )at(A; H )a(^)5t(A; G ) 7 (A; G )) m 

a,f,S 

+v 2 H v 2 g E <« t (^)/?(^)7 t (A:G)7(^)/3 t (^)a(A;H)7 t (A: G )7(A: G )>" 

+ V H V G 

Y,{a\k H )a{k H )5\k G )5{k G )a\k H )a{k H )5\k G )5{k G )) m 
= X + Y 1 + Y 2 + Z . 



(30) 
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Note that we have decomposed each operator into diagonal and off-diagonal parts. We have used 
the condition that the variance of the diagonal matrix elements is twice that of the off-diagonal 
matrix elements in the defining spaces to convert the restricted summations into unrestricted sum- 
mations appropriately to obtain the four terms in the RHS of Eq. (30). Following [47, 49, 51] and 
applying unitary decomposition to 7<5^ (also 5rf) in the first two terms and aft^ (also ftofi) in the 
third term we get X, Y\ and Y 2 . To make things clear, we will discuss the derivation for X term 
in detail before proceeding further. Applying unitary decomposition to the operators ^(k G )S(k G ) 
and ^{k G )8\k G ) using Eq. (29), we have 

k ° / \ 2 

X = v%v 2 a J2 E U-J (»\k H )/3(k H )^ s (s)/3\k H )a(k H )^ s (s)) m . (31) 



Contracting the operators ftft^ across F's using Eq. (21) and operators a^a across T using Eq. 
(20) gives, 

*-44 g (::;) 2 nr) (v) 5 • «> 

Inversion of the equation, 

k G / \ 2 

E <7 t (*G)*(*o)* t (* )7(*o)> m = QH = E ( r - J E (^-M*))™ > (33) 

7,<5 s=0 ^ ' 7,<5 

gives, 



(34) 



For the average required in Eq. (33), we have 

Q(m) = £ (^(k G )5(k G )5\k G Mk G )) m = {^ + k ^) (™) • (35) 

Simplifying Eq. (34) using Eq. (35) and using the result in Eq. (32) along with the series summa- 
tion [47] 

' (_l)*- fl (jy_ t -fc G )!(fc G + t)! fc g !(AT-fc G -g)! /fc G \ /TV + 1\ 
^ (s-f)! (t!) 2 (JV-s-f + 1)! (iV+l-s)! V^/V s /' 1 
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t=o 



the expression for X is, 

X = v 2 H v 2 G F(m, N, k H , k G ) 



' m — s\ 2 ( 'm + kfj — s\ I m — s\ ( m\ ( m\ (N + 1 



^ , j , . x ^ ( m — s\ [m + ha 



N -2s+l [N - s\~ 1 (k G > " 



kn \s \s \ s ) (3V) 



N - s + l \ k G J \s 
Although not obvious, X has kn -H- ka symmetry. This is easy to verify for k H ,kc < 2. In the 
large N limit, Yi, Y 2 and Z are neglected as X will make the dominant contribution; see [51] for 
details on Yi, Y 2 and Z. Thus, in all the applications, we use 

(H(k H )G(k G )H(k H )G(k G )) m = X = v 2 H v% F(m, N, k H , k G ) . (38) 



An immediate application of these averages is in evaluating the fourth order average (H 4 (kH)} m - 
There will be three different correlation patterns that will contribute to this average in the binary 
correlation approximation (we must correlate in pairs the operators for all moments of order > 2), 



{H^k H )) m = ( H(k H )H[ k H )H (k H )H[ k H ))' 



+ { H[k H )H^k H )H^k H )H[ k H )) m m 



+ ( H(k H )H^ H )H(k H )H^ k H )y 



In Eq. (39), we denote the binary correlated pairs of operators with the symbol HH . The first two 
terms on the RHS of Eq. (39) are equal due to cyclic invariance and follow from Eq. (23), 



( H[k H )H[ k H )H [k H )H[ k H )) m = { H{k H )H^k H )H^ H )H {k H )) r 



(40) 



{H\k H )Y 

Similarly, the third term on the RHS of Eq. (39) follows from Eq. (38), 



{H (^k H )H(k H )H[ k H )H(k H )) m = v% F(m, N, k H , k H ) . (41) 



Combining Eqs. (39), (40) and (41), {H A (k H )) m is given by, 

i 2 



{H\k H )) m = v A H [2 {T(m, N, k H )} + F(m, N, k H , k H )] . (42) 
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B. Results for two unitary orbits 



In the NDBD applications (also (5 decay), fourth order traces over two orbit configurations 
are needed. Let us consider m particles in two orbits with number of sp states being Ni and 
N 2 respectively. Now the m -particle space can be divided into configurations (m 1; m 2 ) with mi 
particles in the #1 orbit and m 2 particles in the #2 orbit such that m = mi + m 2 . Considering the 
operator H with fixed body ranks in mi and m 2 spaces such that (mi, m 2 ) are preserved by this 
operators, the general form for H is, 

H{k H )= [ v H l5 {h3)} a\(i)Mi)4UMj) ■ (43) 

Now, it is seen that, in the dilute limit, 



2(U„\\ mi ' m2 



mi ,7712 



= E E («I(OA(07lO> 2 (j)^««i«40')72(i) 



E E(^w^«^«^«) E(^o> 2 (j)4(j)7 2 (j) 

i+j=kn otfi 7,<5 

E T{m u N u i) T(m 2 ,N 2 ,j) . 

i+j=kn 



rr/2 



(44) 



Note that v 2 H (i,j) = [v^ l5 {i,j)] 2 and T's are defined by Eqs. (24) and (25). The ensem- 
ble is defined such that v^ l5 {i,j) are independent Gaussian random variables with zero center 
and the variances depend only on the indices i and j. Similarly, with two operators H and G 
(with body ranks k H and k G respectively) that are independent and both preserving (m 1: m 2 ), 



(H(k H )G(k G )H(k H )G(k G )} mi ' m * is given by, 



(H(k H )G(k G )H(k H )G(k G )) m ^ = 

(45) 

E v H(hj) v g(^ m ) F^m^N^^t) F(m 2 ,N 2 ,j,u) , 

i+j=kn, t+u=ka 

15 



and therefore, 

2 



(H\k H )) mi ' m2 = 2 



v 2 H {i,j)T{m u N u i)T{m 2 ,N 2 ,j) 



(46) 



+ Yl v H(hj) v H(t, u ) F(mi,JVi,i,*) F(m 2 ,N 2 ,j,u) . 

i+j=kH, t+u=k H 



Now we apply the formulation given here to derive the formulas for the second and fourth order 
cumulants denning p Uv _g.Q.H • 



V. BINARY CORRELATION RESULTS FOR THE BIVARIATE CORRELATION COEFFI- 
CIENT AND FOURTH ORDER CUMULANTS FOR THE TRANSITION STRENGTH DENSITY 
FOR NDBD 

A. Transition matrix elements and bivariate strength densities 

Our purpose here is to establish that for the 0v(3~(3~ decay (also for (3 decay), transition strength 
densities, locally, are close to bivariate Gaussian form and also to derive a formula for the corre- 
sponding bivariate correlation coefficient. With space #1 denoting protons and similarly space #2 
neutrons, the general form of the transition operator O is, 

0{k ) = J2 v o( k o) ll(ko)S 2 (k ) ; k a = 2 for NDBD . (47) 

Therefore, in order to derive the form for the transition strength densities generated by O, it is 
necessary to deal with two-orbit configurations denoted by (mi, m 2 ), where mi is the number of 
particles in the first orbit (protons) and m 2 in the second orbit (neutrons). Now, the transition 
strength density Io(Ei, Ef) is 

l (Ei,E f ) 

(48) 

= /Kw)/^) \(( mi ,m 2 ) f E f | O | {m^m^E^ 2 I^^iEi) , 
and the corresponding bivariate moments are 



M PQ (( mi ,m 2 ) t ) = (0^k )HQ(k H )0(k )HP(k H ))^ m ^ . (49) 
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Note that M are in general non-central and non-normalized moments. The general form of 
the operator H{k H ) is given by Eq. (43) and it preserves (mi,m 2 )i's. However, O and its 
hermitian conjugate do not preserve (m 1: m 2 ) i.e., 0(k ) \m 1: m 2 ) = \m 1 + k ,m 2 — k a ) 
and 0\ko) |m 1 ,m 2 ) = \m 1 — k ,m 2 + k ). Thus, given a (m 1 ,m 2 )j for an initial state, the 
(mi,m 2 )/ for the final state generated by the action of O is uniquely defined and therefore, the 
bivariate moments defined by Eq. (49) are proper bivariate moments and they are defined by 
the initial (mi, m 2 )j. For completeness, let us mention that given the marginal centroids (e, h e/), 
widths (tJi, 07) and the bivariate correlation coefficient 0>™, the normalized bivariate Gaussian is 
defined by, 

Pbiw-g-,o( E ^ E f) = Pbiy-g-,o( E ^ E f e *' e /' a *> °7> Cw«) 
1 

~ WFU (50) 



B. Formulas for the bivariate moments 

Using binary correlation approximation, we derive formulas for the first four moments M P q ( (m x 
, m 2 )i),P + Q < 4 of I ^ m ^ m ^!( K E i , E f ) for any k a by representing H(k H ) and 0(k ) op- 
erators by independent EGOEs and assuming H{k H ) is a /c^-body operator preserving (m 1 ,m 2 )'s. 
Note that the ensemble averaged /c^-particle matrix elements of H{k H ) are v 2 H (i,j) with i + j > = 
k H [see Eq. (43)] and similarly the ensemble average of (v^) 2 is v 2 . From now on, we use 
(mi, m 2 )i = (mi, m 2 ). Using Eq. (47) and applying the basic rules given by Eqs. (18) and (19), 
we have 



M 00 (mi,m 2 ) = (0^k )0(k )} 



mi,m2 



E (4(^)7i(A:o)7i t (A:o)5 2 (A;a)) mi ' m2 (51) 

7,5 



mA /m 2 
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Trivially, M 10 (m 1 ,m 2 ) and M 01 (m 1 ,m 2 ) will be zero as H{k H ) is represented by EGOE (/;#). 
Thus, Mpg(mi,m 2 ) are central moments. Moreover, by definition, all the odd-order moments, 
i.e., MpQ(mi,m 2 ) with mod (P + Q,2) ^ 0, will be zero. Now, the M u is given by, 



M 11 (m 1 ,m 2 ) = (0^k )H(k H )0(k )H(k H )) 



mi,m2 



mi 



u& ^ «lr(i,j) {ll(k )a 1 (i)f3\(i) ll (k )(ii(i)a\{i) t 

01,01,012,02,71,82; i+j=k H 

(52) 

x (5 2 (ko)a 2 (j)Pt(j)5t(ko)P2(j)al(j)" 



Then, contracting over the 7*7 and 55^ operators, respectively in the first and second traces in Eq. 
(52) using Eqs. (20) and (21) appropriately, we have 

M 11 (m 1 ,m 2 ) = vl £ (VOCV') 

(53) 

x T(m 1 ,N 1 ,i)T(m 2 ,N 2 ,j) . 

Note that the formulas for the functions T(- • • )'s appearing in Eq. (53) are given by Eqs. (24), 
(25) and (26). Similarly, the functions F(- • • )'s appearing ahead are given by Eqs. (28) and (37). 
For the marginal variances, we have 



M 20 (m 1 ,m 2 ) = (0^k )0(ko)H 2 (k H )) 



mi,m.2 



= M 00 ( mi ,m 2 ) (H 2 (k H )) 



mi,m,2 



(54) 



M 02 ( mi ,m 2 ) = (Oi(k )H 2 (k H )0(k )} 



mi ,7712 



= M 00 ( mi ,m 2 ) (H\k H )) 



m 1 +k ,m 2 -k 



In Eq. (54), the ensemble averages of /J 2 (/c#) are given by Eq. (44). Now, the correlation 
coefficient (, biv is 

r , \ Mu (7711,7712) , <c . 

C^(m 1 ,m 2 )= = (55) 

\J M 20 (m 1 ,m 2 ) M 02 (m 1 ,m 2 ) 

Clearly, £ Wt; will be independent of v 2 Q . 
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Proceeding further, we derive formulas for the fourth order moments M PQ , P + Q = 4. The 
results are as follows. Firstly, for (PQ) = (40) and (04), we have 



M 40 (m 1 ,m 2 ) = (01(k )0(ko)H4(k H )) 



mi ,7712 



M o(mi,m 2 ) (H\k H )Y 



(56) 



M 04 (m 1 ,m 2 ) = (0^k )H\k H )0{k )) 



mi ,7712 



= M 00 (m l! m 2 ) (tf^^p+ta-fco . 
In Eq. (56), the ensemble averages of H 4 (k H ) are given by Eq. (46). For (PQ) = (31), we have 



M 31 (m 1 ,m 2 ) = (0t(A; o )#(A;H)0(M# 3 (A;H)> 



mi, 7772 



(0^ko) H{k H )0{ko)H^ H ) H{k H )H {k H )Y 



+ {0^ko) H^k H )0{ko)H^k H )H{k H )H^ k H )y 



+ (0^(k )H(k H )0(k )H(k H )H(k H )H(k H )) n 



(57) 



j 



First and last terms on RHS of Eq. (57) are simple as HH can be taken out of the average and 
then we are left with a term similar to M n (mi ,m 2 ). For the second term, the and O operators 
are contracted across H operator using Eqs. (20) and (21) and then one is left with an average of 
the form (HGHG). These will give the final formula, 



M 3l (m u m 2 ) = 2 (H^(k H )) mi ' m2 M 11 (m 1 ,m 2 ) 



+ (OKko)H ^k H )0(k )H(,k H )H[k H )H^ k H )) m ^ 



(58) 

2 (H*(k H )) m ^ M u (m u m 2 )+v 2 £ v 2 H (i,j) v 2 H (t,u) 

i+j=k H ,t+u=k H 



y[ m \ o ') r\ o l ) F{n, 1 .X x .i.t)F{m 2 .X 2 .j.u) 
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Similarly, we have 



mi,m2 



M 13 (m 1 ,m 2 ) = (Crt(k )H*(k H )0(ko)H(k H )) 



(0^k )H ^ H )H{ k H ) H{k H )0{ko)H {k H )y 



+ (OHko)H ^ H )H^ H )H^H)0(ko)H^ k H )) mi ^ 



+ (Ql(ko)H [k H )H(k H )H(k H )0(k )H[ k H )} 



rai,m2 



= 2 (H 2 (k H )} mi+ko ' m2 - ko M n (m 1 ,m 2 ) 



+v 2 o v 2 H (i,j)v 2 H (t,u) G(t,u) 

i+j=k H ,t+u=k H 



X 



(59) 



m i — ko — t + i\ ( m 1 + k a — t\ f m 2 — u + k + j\ f m 2 — k — u 
i ) V i ) V 3 J \ 3 



G(t,u) = U ^T(mi,iVi,t) T(m 2 ,iV 2!M ) . 

In Eq. (59), the first and last terms can be evaluated by first calculating the H 2 average over the 
intermediate states \m\ + k , m 2 — ko) and then the remaining part is similar to Mn(mi, m 2 ). 
Also, the second average is evaluated by first contracting the two correlated H's that are between 
and O operators (see the contraction symbol for clarity) and then one is again left with a term 
similar to M 11 (m 1 , m 2 ). Finally, M 22 (toi, m 2 ) is given by, 



M 22 (m 1 ,m 2 ) = (0^k )H\k H )0{k )W{k H )) 



mi,m2 



= {0^ko)H ^k H )H {k H )0{ko) H{k H )H {k H )) 



+ {OKko) H{k H )H{k H )0{k )H {k H )H{k H )) 



+ (OKka) H(k H )H(k H )0(k )H(k H )H^ k H )) 
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= M 00 (m 1 ,m 2 ) ( H i(k H )) mi+ko ' m2 - ko (H 2 (k H )) mi ' m2 

i+j=k H ,t+U=k H V V / \ / 

(60) 

x [FimuN^t) F(m 2 ,N 2 ,j,u) 

+T(m 1; N u i) T(m u N u t) T(m 2 , N 2 ,j) T(m 2 , N 2 , u)} . 
In Eq. (60), the first term is evaluated by first calculating the H 2 average (for the H 2 between 
C?t and O operators) over the intermediate state \rn-y + k , m 2 — ko) and then one is left with 
product of averages of H 2 and O^O operators. For the third term, first the and O operators 
are contracted across H 2 operator and then we are left with average of the form (H 2 ) x (H 2 ). 
Similarly, for the second term, after contracting the and O operators across H 2 operator, we 
are left with an average of the form (HGHG) . 



C. Numerical results for bivariate correlation coefficient and fourth order cumulants 

Firstly, given the M P Q(mi, m 2 ), the normalized central moments M P q are M P q = M P q/M 00 . 
Then, the scaled moments M P q are 

m pq = MpQ p/ ( 2 mi,m2) or^ p + Q>2. (61) 

Now the fourth order cumulants are [52], 

fc 4 o(mi, m 2 ) = M 40 (mi, m 2 ) - 3 , k Gi (mi,m 2 ) = M 04 (mi, m 2 ) - 3 , 

fc 3 i(mi,m 2 ) = M 31 (m 1 ,m 2 ) - 3 M n (m 1 ,m 2 ) , 

^ ^ (62) 

k lz (m 1 ,m 2 ) = M n (m l ,m 2 ) - 3 M ll (m 1 ,m 2 ) , 

k 22 (m 1 ,m 2 ) = M 22 (m l ,m 2 ) - 2 M 1 2 1 (mi,m 2 ) - 1 . 

Assuming vjj(i, j) defining H{2) are independent of (i, j) so that ( biv is independent of v 2 H , we 
have calculated the value of C,uv with ko = 2 for several 0u(3~ (3~ decay candidate nuclei using Eq. 
(55) along with Eqs. (51), (53), (54) and (44). For the function T(- • • ), we use Eq. (24). Note that 
vjj(i, j) correspond to the variance of two-particle matrix elements from the p — p (i = 2, j = 0), 
n — n (i = 0, j = 2) and p — n (i — 1, j — 1) interactions. Results are given in Table I. It is seen 
that ( biv ~ 0.6 — 0.8. It is important to mention that = for GOE. Therefore, the transition 
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TABLE I. Correlation coefficients (biv{ m i> m 2) f° r some nuclei with k@ = 2 as appropriate for Oz//3~/3~ 
decay operator. Note that space #1 is for protons and space #2 for neutrons. The configuration spaces 
corresponding to N\ or N 2 = 20, 22, 30, 32, 44 and 58 are rsf, r%g, r±g, r^h, r$i, and r^j, respectively 
with / = 1 f 7/2 , g = l g 9 / 2 , h = 1 h ll/2 , i = 1 h?,/2, 3 = 1 Jib/2, n = 1 h/2 2 P3/2 2 Pi/2, n = 1 g 7/2 2 d 5/2 2 d 3/2 
3 s 1/2 , r 5 = 1 h g/2 2 f 7 /2 2 f5/2 3 P3/2 3 Pi/2 and r 6 = 1 i 11/2 2 g 9/2 2 g 7/2 3 d 5/2 3 d 3/2 4 s 1/2 . See text for details. 



Nuclei 




m\ 


N 2 


m 2 


(biv(mi,m 2 ) 


32 Ge 44 


22 


4 


22 


16 


0.64 


34Se48 


22 


6 


22 


20 


0.6 


JS°M068 


30 


2 


32 


8 


0.57 


i 8 Te 76 


32 


2 


32 


26 


0.62 


5 fTe 78 


32 


2 


32 


28 


0.58 


60°Nd 90 


32 


10 


44 


8 


0.72 


g2 4 Sm 9 2 


32 


12 


44 


10 


0.76 


JfWioe 


32 


24 


44 


24 


0.77 


ii 8 ui46 


44 


10 


58 


20 


0.83 



strength density will be narrow in (Ei, Ef) plane. In order to establish the bivariate Gaussian form 
for the 0u(3~ (3~ decay transition strength density, we have examined kpQ, P + Q = 4. For a good 
bivariate Gaussian, \k PQ \ < 0.3. Using Eqs. (51), (53), (54), (56), (58)-(62) along with Eqs. (44) 
and (46), we have calculated the cumulants kpQ(mi : m 2 ), P + Q = 4. These involve T(- • • ) and 
F(- ■ ■ ) functions. For set #1 calculations in Table II, we use Eq. (24) for T(- • • ) and Eq. (37) for 
F(- ■ ■). For the set #2 calculations, shown in 'brackets' in Table II, weuseEq. (25) for T(- ■ ■ ), Eq. 
(28) for F(- ■ ■ ) and replace everywhere ( mi s +r ) — > { Ni ) for any (r, s) with i — 1,2. Then we have 
the strict dilute limit. We show in Table II, bivariate cumulants for five heavy nuclei for both sets 
of calculations and they clearly establish that bivariate Gaussian is a good approximation (similar 
tests are made for j3 decay operator in Appendix A). We have also examined this analytically in 
the dilute limit with Ni,N 2 — > oo and assuming vjj(i,j) independent of With these, we 

have expanded k PQ in powers of l/m 1 and l/m 2 using Mathematica. It is seen that all the k PQ , 
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TABLE II. Cumulants kpQ, P + Q = 4 for some nuclei listed in Table I. The numbers in the brackets 
are for the strict dilute limit as explained in the text. Just as in the construction of Table I, we use vfj(i,j) 
independent of (i, j). See Table I and text for details. 



Nuclei Ni mi N 2 m 2 k± fe 4 k 13 k 31 k 22 



%°Mo 58 30 


2 


32 


8 


-0.45(- 


-0.39) 


-0.42(- 


-0.38) 


-0.24(- 


-0.23) 


-0.26(- 


-0.25) 


-0.20(- 


-0.22) 


i[> Nd 9 o 32 


10 


44 


8 


-0.27(- 


-0.22) 


-0.29(- 


-0.23) 


-0.22(- 


-0.18) 


-0.20(- 


-0.17) 


-0.19(- 


-0.18) 


g2 4 Smg2 32 


12 


44 


10 


-0.24(- 


-0.18) 


-0.25(- 


-0.18) 


-0.19(- 


-0.15) 


-0.18(- 


-0.15) 


-0.17(- 


-0.15) 


74°Wi06 32 


24 


44 


24 


-0.19(- 


-0.08) 


-0.20(- 


-0.08) 


-0.17(- 


-0.08) 


-0.15(- 


-0.08) 


-0.15(- 


-0.08) 


1| 8 U 146 44 


10 


58 


20 


-0.18(- 


-0.13) 


-0.18(- 


-0.13) 


-0.15(- 


-0.11) 


-0.15(- 


-0.11) 


-0.13(- 


-0.11) 



P + Q = 4 behave as, 

kpQ = — + o(±)+o(?£)+... . (63) 
mi \ m iJ \ m iJ 

3 /2 

Therefore, for mi >> 1 and m 2 « m/ , the strength density approaches bivariate Gaussian 
form in general. It is important to recall that the strong dependence on rxi\ in Eq. (63) is due to 
the nature of the operator O i.e., 0(k ) \m 1: m 2 ) = \mi + k , rn 2 — ho). Thus, we conclude that 
bivariate Gaussian form is a good approximation for 0u/3~/3~ decay transition strength densities. 
With this, one can apply the formulation given in Sec. Ill with the bivariate correlation coefficient 
( biv given by Eqs. (55), (54) and (53). The values given by the two-orbit binary correlation theory 
for (bi V can be used as starting values in practical calculations. 

VI. CONCLUSIONS 

To summarize, by extending the binary correlation approximation method for spinless embed- 
ded A:-body ensembles to ensembles with proton-neutron degrees of freedom that involves traces 
involving product of powers of two different operators over two-orbit configurations (either the 
operators preserve the two-orbit symmetry or change a two-orbit configuration to a unique final 
configuration), we have established that the transition strength density generated by the two-body 
part of the Hamiltonian is a bivariate Gaussian for transition operators 0{ko) that change k 
number of neutrons to k G number of protons. Towards this end, we have derived formulas for 
the fourth order cumulants of the transition strength density and calculated their values for some 
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TABLE III. Correlation coefficients Cftii;( m i> m 2) and cumulants kpQ, P + Q = 4 for some nuclei relevant for @ decay. The first four nuclei in the table 
are relevant for (3~ transitions, next four nuclei are relevant for electron capture and the last two nuclei are relevant for (3 + transitions. The numbers 
in the brackets for kpQ are for the strict dilute limit as in Table II. We assume v 2 H {i,j) are independent of as used in the calculations generating 
Tables I and II. Here, m\ = m p , m 2 = m n for the first four nuclei and mi = m n , m 2 = m v for the next six nuclei. See text for details. 



Nuclei 


Ni 


mi 


N 2 


m 2 


(bi v (mi,m 2 ) 


ho 


k 


04 


k 


13 


k 


31 


k 


22 


27C035 


20 


7 


30 


15 


0.72 


-0.26(-0.18) 


-0.27( 


-0.18) 


-0.24( 


-0.16) 


-0.23( 


-0.16) 


-0.22( 


-0.16) 


27C037 


20 


7 


30 


17 


0.73 


-0.27(-0.16) 


-0.27 ( 


-0.16) 


-0.24( 


-0.15) 


-0.23( 


-0.15) 


-0.21( 


-0.15) 


26 Fe 36 


20 


6 


30 


16 


0.72 


-0.28(-0.18) 


-0.28( 


-0.18) 


-0.24( 


-0.16) 


-0.24( 


-0.16) 


-0.22( 


-0.16) 


lNi 40 


20 


8 


30 


20 


0.72 


-0.27(-0.14) 


-0.27( 


-0.14) 


-0.24( 


-0.13) 


-0.23( 


-0.13) 


-0.21( 


-0.13) 


llGe 33 


36 


5 


36 


4 


0.55 


-0.45(-0.41) 


-0.46( 


-0.42) 


-0.35( 


-0.33) 


-0.34( 


-0.32) 


-0.34( 


-0.34) 


ilSe 35 


36 


7 


36 


6 


0.66 


-0.36(-0.29) 


-0.34( 


-0.30) 


-0.28( 


-0.25) 


-0.28( 


-0.25) 


-0.27 ( 


-0.25) 


3|Kr 37 


36 


9 


36 


8 


0.72 


-0.28(-0.23) 


-0.28( 


-0.23) 


-0.24( 


-0.20) 


-0.24( 


-0.20) 


-0.23( 


-0.20) 


38Sr 39 


36 


11 


36 


10 


0.76 


-0.24(-0.19) 


-0.24( 


-0.19) 


-0.21( 


-0.17) 


-0.21( 


-0.17) 


-0.20( 


-0.17) 


iMo 43 


36 


15 


36 


14 


0.79 


-0.20(-0.14) 


-0.21( 


-0.14) 


-0.19( 


-0.13) 


-0.18( 


-0.13) 


-0.17( 


-0.13) 


iPd47 


36 


19 


36 


18 


0.80 


-0.19(-0.11) 


-0.19( 


-0.11) 


-0.18( 


-0.10) 


-0.17( 


-0.10) 


-0.16( 


-0.10) 



realistic examples; they are found to vary from ~ —0.4 to —0.1. It is important to mention that 
the embedding algebra for the EGOEs used is U (N) D U (N p ) © U (N n ) [p denotes 'protons' and 
n denotes 'neutrons'] with the Hamiltonian preserving the symmetry and the transition operator 
breaking the symmetry in a particular way. We have also derived a formula for the fourth order 
trace defining the correlation coefficient of the bivariate transition strength density for the transi- 
tion operator relevant for 0isj3~l3~ decay. For nuclei from 76 Ge to 238 U, the bivariate correlation 
coefficient is found to vary from ~ 0.6 — 0.8 and these values can be used as a starting point for 
calculating nuclear transition matrix elements for NDBD using the spectral distribution method 
outlined in Sec. III. In future, it is important to test the approximations leading to Eq. (14) using 
shell model examples. Although spectral distribution method is expected to be valid in the chaotic 
domain of the spectrum (usually away from the ground state), it remains to be tested how well the 
method applies to the calculation of NTME for NDBD. In the past, the theory has been applied 
successfully for occupancies near the ground state [19, 53-55] and also it is shown that in the level 
density analysis of heavy nuclei [44] that the theory extends close to the ground state. In the near 
future, applications will be carried out for NTME for some heavy nuclei ( 100 Mo, 154 Sm, 150 Nd, 
186 W, 238 U). 
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APPENDIX A 

For completeness, we have also calculated the correlation coefficient and fourth order moments 
for the transition operator relevant for j3 decay [k a = 1 in Eq. (47)]. Results are given in Table 
III. For the first four nuclei (they are (3~ candidates) in the table, N = Z = 20 is used as the core. 
Here, TVi corresponds to 1 f 7/2 1 f 5/2 2 P3/2 2 Pi/2 and iV 2 corresponds to 1 f 7/2 Vs/2 2 Ps/2 2 Pi/2 1 9%/2- 
Similarly, for the remaining six nuclei (they are electron capture or (3 + candidates), N = Z = 28 
and Ni and N 2 correspond to 1 / 5 / 2 2 P3/2 2 Pi/2 1 Q9/2 1 9i/2 2 d$/ 2 . The fourth order cumulants values 
presented in Table III confirm that the bivariate Gaussian form is a good approximation for j3 decay 
transition strength densities. Results in Table III justify the assumptions made in [39, 56] where 



25 



spectral distribution method is applied, with the correlation coefficients in the correct range, to 
calculate the ft decay rates for nuclei relevant for pre-supernovae evolution. 
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